Quasi-maximum likelihood estimation of periodic GARCH processes 
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Abstract 

This paper establishes the strong consistency and asymptotic normality of the quasi-maximum like- 
lihood estimator (QMLE) for a GARCH process with periodically time-varying parameters. We first 
give a necessary and sufficient condition for the existence of a strictly periodically stationary solution for 
the periodic GARCH (P-GARCH) equation. As a result, it is shown that the moment of some positive 
order of the P-GARCH solution is finite, under which we prove the strong consistency and asymptotic 
normality (CAN) of the QMLE without any condition on the moments of the underlying process. 
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Introduction 

5-H 



Periodic GARCH (P-GARCH) processes introduced by Bollerslev and Ghysels (1996) have proved useful 
and appropriate for modeling many time series encountered in practice, which are characterized by a 
stochastic conditional variance with periodic dynamics (see also Franses and Paap 2000; 2004; Ghysels and 
Osborn, 2001 and the references therein). As for the standard GARCH model, the P-GARCH process may 
be seen as a non Gaussian white noise but whose conditional variance follows a linear periodic ARM A 
(PARMA) dynamics in terms of the squared process, implying a periodic structure for the underlying 
process itself. In contrast with periodic ARM A models which may be written as an equivalent vector 
ARM A form (Tiao and Grupe, 1980), there is no direct correspondence between the P-GARCH model 
and the multivariate GARCH one. In fact, any P-GARCH process can be written as only a weak 
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multivariate GARCH model (see Drost and Nijman (1993) for the definition of weak GARCH) meaning 
that the study of the P-GARCH may not be trivially deduced from the existing multivariate GARCH 
theory and thus it constitutes a useful and interesting task. 

Since their introduction, P-GARCH models have been fairly considered in the literature. Bollerslev and 
Ghysels (1996) have studied the second-order periodic stationarity of the P-GARCH model and derived 
the corresponding QMLE but without studying its asymptotic properties. They have successfully applied 
the P-GARCH modeling for certain real data. Franses and Paap (2000) have considered the P-GARCH 
process as an innovation of a more general periodic ARM A model with P-GARCH error. Important 
applications of such processes may be found in Ghysels and Osborn (2001) and Franses and Paap (2004). 
On the other hand, some probabilistic properties such as strict periodic stationarity, existence of higher 
order moment and geometric ergodicity have been studied recently by Bibi and Aknouche (2006). Apart 
from the mentioned works, it seems that there is no result concerning asymptotic inference about such 
models. However this problem has been exhaustively studied recently for the standard GARCH case. 
Indeed, a considerable amount of research has been executed for studying the asymptotic properties of the 
QMLE for GARCH processes. This research, including contributions by authors such as Lee and Hansen 
(1994), Lumsdaine (1996), Boussama (2000), Berkes et al. (2003), Ling and McAleer (2003), Berkes and 
Horvath (2004), and Francq and Zakoian (2004), is aimed at establishing consistency and asymptotic 
normality of the QMLE for GARCH processes with weak conditions on the moment as this makes an 
undesirable restriction on the parameter space. Francq and Zakoian (2004) 's result seems to be the weaker 
one obtained at present. 

In the spirit of the latter work, the main goal of this paper is to establish the strong consistency 
and asymptotic normality of the QMLE for the P-GARCH process. We first show that the sufficient 
condition for strict periodic stationarity of the P-GARCH process given by Bibi and Aknouche (2006) 
is also necessary. As a consequence, it will be shown that the moment of some positive order of the 
periodic stationary solution is finite. This result will be exploited in order to obtain strong consistency and 
asymptotic normality for the QMLE irrespective of any moment requirements for the underlying process. 

The rest of this paper is organized as follows. Section 2 gives a necessary and sufficient condition for 
the strict periodic stationarity for the P-GARCH process and some important consequences thereof. The 
main results concerning strong consistency and asymptotic normality of the QMLE are found in section 
3 while the proofs of the main results are left in the appendix. 

Through the rest of this paper we need the use of the following notations: the symbols " and 
"a.s. stand for convergence in law and almost sure convergence respectively. A matrix A of order m x n 
is denoted by A mxn and the k x k identity matrix by Ikxk- p{A) denotes the spectral radius, i.e. the 
maximum modulus of the eigenvalues of squared matrix A, A' stands for the transpose of A and the 
vectorial relation A> B means that each element of A is greater than the corresponding element of B. 
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2 Strict periodic stationarity of the P-GARCH and some of its impor- 
tant consequences 



This paper studies the periodic GARCH process (yt,t £ Z) of orders p and q and period S > 1, that is a 
solution to the stochastic difference equation 



Vt = VhtVt 

u , <f 2 , £ o z, . *€Z = {0,±1,±2,...}, 

i=i j=i 



(2.1) 



where G Z) is a sequence of independent and identically distributed (i.i.d.) random variables defined 
on a probability space (fi, „4, P) such that -B (%) = 0, E = 1 and rj k is independent of y t for all fe > t. 
The parameters w t , a t: i, and /3 t ■ are periodic in t with period 5 (i.e., 0Jt+kS = w t) a t+ks,i = a t,i and 
Pt+ks,j = Pt,j) such that uj t > 0, ot ti i > and f3 tj - > 0, « = 1, ...,<?, J = 1, —,p, for all k,t £ Z, so that by 
setting t = v + Sn, model (2.1) may be equivalently written as 



Vv+Sn = y/K+SnVv+Sn 

h v+Sn =UJ V + a V ,iy v+Sn -i + E Pv,j h v+Sn-j 
i=l j=l 



(2.2) 



highlighting thus the periodicity in the model. In the difference equation (2.2) y v +sn may refer to (yt) 
during the v — th 'season' of the year n and u v , a V; ±, together with P v i, P v p are the model 
coefficients at season f . The orders p and g can also be considered periodic in time, say pt and qt in 
order to reduce the number of parameters to be estimated. However, in this case a constrained estimation 
procedure, which is not the objective of this paper, would be needed. 

A standard approach for studying the structure of an ARCH-type model is to write it as a random 
coefficient representation. Indeed, it is easy to write model (2.1) in term of the squared process as follows 



yt = utm + E ®t,iVtVt~i + E PtjVtht-j 



i=i j=i 

q p 

h t = uj t +J2 ot t ,iVt-i + E Ptjht-j, 



t G Z, 



(2.3) 



i=l j=l 

which is ready to be cast in a stochastic recurrence equation with random coefficients. 
Indeed, defining the (p + ^-random vectors Y t = (y% , h t , ht- p +i)' and 

B t = (wtrjt, 0( g _i) x i) ^ti O(p-i)xi) ' ("mxn stands for the null matrix of order mxn) together with the 
(p + q) x (p + q) random matrix A t given by 



/ 



At 



(Xt,lVt a t,2Vt ■ ■ ■ at,q-lVt ®t,qVt 



0t,iVt 0t,2Vt ■ ■ ■ Pt, P -iVt Pt, P Vt \ 



I (q-l)x(q-l) 
&t,l OL t ,2 ■ ■ ■ <*t,q-l 
%-l)xq 







(9-1) Xl 

a t , q 







(9-1) xp 



Pt,l P 



t,2 



Pt,p-1 P 



t,p 



J (p-l)x(p-l) %-l)xl / 
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one can rewrite model (2.3) in the following generalized AR model 

Y t = AM-! +Bt,t€ Z, (2.4) 

which differs from the standard formulation studied by Bougerol and Picard (1992a, 19926) in that the co- 
efficients (A t ,B t ) are rather independent and periodically distributed (i.p.d.). Moreover, their expectations 
E {A t } and E {B t } defined element-wise are 5-periodic in time. 

Now, we focus on the stochastic recurrence equation (2.4) and the existence of a strictly periodically 
stationary (henceforth s.p.s.) solution thereof. Recall that the process (Y t , t G Z) is said to be s.p.s. (with 
period S > 1) if the distribution of (Y tl , Y t2 , Y tn )' is the same as that of (Y tl+ sh, Y t2+S h, — , Y tn+S h)' for 
all n > 1 and h, £i,t2, ■■■,t n G Z. Furthermore, it is called periodically ergodic (cycloergodic, cf, Boyles and 
Gardner, 1983) if for all Borel set B and all integer m 

1 n 

- Ifi {Yy+st, Yy+i+St, Y v+m +st) — ► P({Y v ,Y v+ i, ...,Y v+m ) G B) , a.s. as n -»■ oo, 
n t=i 

for all 1 < v < S, where denotes the indicator function of the set B. It follows immediately that a 
standard ergodic process is periodically ergodic with period S = 1. As for the stationary case (Billingsley, 
1996 Theorem 36.4), periodic ergodicity is closed under certain transformations. In particular if (et, t G Z) 
is s.p.s. and periodically ergodic and if (Yt, t G Z) is given by Yj = / (..., £t,et+i, •••) where / is a measurable 
function from R°° to R, then (Yj, i G Z) is also s.p.s. and periodically ergodic (Aknouche and Guerbyenne, 
2007). A periodic analog of the ergodic theorem for stationary sequences can be stated as follows. If 
(Yt, t G Z) is s.p.s. and periodically ergodic and if / is a measurable function from R°° to R such that 
E{f(...,Y t - 1 ,Y t ,Y t+1 ,...)} < oo then for all 1 < v < S, 

1 n 

~y^f {-, Y v+S(t-i), Y v+St,Y v+S (t+i),-) E {f {-, Y v+S(t-i), Y v+St, Y v+S(t+i),-)} i a - s - as n > oo. 
n t=i 

It is well known that the periodic process theory can be translated immediately into the stationarity process 
one by an appropriate transformation. Since the seminal paper by Gladyshev (1961), we know that a 
periodically stationary process (Y t , t G Z) is equivalent to a vector- valued stationary process (Y n , n G Z) 
where 

Hn = Y 2+n 5, Y S+ ns)' , 

meaning that any property about periodic processes translates at once into a corresponding result about 
stationary processes. In our case, the process (Y n ,n G Z) satisfies nonuniquely the stochastic difference 
equation 

Y n = A n Y n _ 1+ B n , nGZ, (2.5) 
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where A n and B_ n are defined by blocks as 



f rX r 
Or xr 



A, 



r 



o r 



A\+nS 
M+nsM+nS 



5-1 

11 A nS+S-j 
v=0 



, B, 



B\+nS 
A2+nsB\ +n S + B 2 +nS 



s rs-k-i 

Yl { 11 A S^v+nS } B k+nS 
fc=l L v=0 



' rSxrS x fc=l ^ v =v J ' rSxl 

with r = p + q and where, as usual, empty products are set equal to one. Therefore, the solution process 
(2.4) is s.p.s. if and only if the corresponding solution of (2.5) is strictly stationary. Similarly, the process 
given by (2.4) is periodically ergodic if and only if the corresponding process solution of (2.5) is ergodic. 
So, we can substitute the study of the properties of model (2.4) to those of model (2.5) if they are much 
easy to obtain. In the sequel we only use model (2.4) as it seems to be simpler. 

As for strict stationarity, the primordial tool in studying strict periodic stationarity is the top Lyapunov 
exponent for i.p.d. random matrices, properties of which can be found in Aknouche and Guerbyenne (2007). 
Let ||.|| be an arbitrary operator norm in M r , the space of real square matrices of dimension r. Then, the 
top Lyapunov exponent associated with the i.p.d. sequence of matrices A := (At, t G Z) is defined when 
Eti^(log + P,||) <oo,by 



r S (A)= inf -E{\og\\A nS A nS _ 1 ...A 1 \ 

n€N* n 



(2.6) 



where log + (x) = max(log (x) , 0). For S = 1, (2.6) reduces to the definition of the top Lyapunov exponent 
for i.i.d. matrices (e.g. Bougerol et Picard, 1992a). It is clear that 7 5 (.) inherits the properties of the 
standard top Lyapunov exponent. In particular, the following inequality 7 5 (A) < J2v=i ^ ll^fll) 
holds with equality for the scalar case r = 1. Furthermore, when the (A t , t G Z) is a sequence of ^-periodic 
nonrandom matrices, then from (2.6) 



7 S (A) = log lim 



lu 

\v=0 



S-v 



1/n 



def 

= logp 



\v=l 



On the other hand, by the i.p.d. property of (A t , t G Z), the sequence ^J^Zq 1 A n s- V , n £ Z^ is i.i.d. to 
which applying the subadditive ergodic theorem (cf., Kingman 1973) or the result of Furstemberg and 
Kesten (1960) it follows that a.s. 



7 5 (A)= lim -logP n5 A n5 _i...^i| 

n^oo n 



(2.7) 



As noted below, properties of model (2.4) can be obtained from those of model (2.5) which is a generalized 
AR with nonnegative i.i.d. coefficients for which a fairly exhaustive theory exists (see Bougerol and Picard 
1992a, 19926 and the references therein). Since {rf t ,t 6 Z) is i.i.d. with finite variance and the parameter 
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are S-periodic, the sequence ((A t ,B t ) ,t€Z) is s.p.s. and periodically ergodic and then the random pair 
{(A t ,B_ t ) ,t G Z) is stationary and ergodic. Moreover, because Eli ^ 0°g + ||A>||) < Eli ^ (II A; II) an d 



Eli # (log + II^H) < Eli # (II^H), then Eli # 0og + P,||) and Eli # (log + ||A,||) together with 
E (log + HAqII) and £ (log + HBqH) are finite. Thus from Theorem 2.3 of Bougerol and Picard (19926) 



a necessary and sufficient condition for model (2.5) to possess a nonanticipative (future-independent) 
strictly stationary solution is that the top Lyapunov exponent associated with A = (A t ,t £ Z), that is 
7 (A) := infteN* jE {(log || ; A t .A t _ 1 ... ; A 1 |h } is strictly negative. 

It must be noted that the latter result is not suitable since it gives a condition about the transformed 
non-periodic model (2.5), not the original periodic one (2.4), which is the objective of this paper. For- 
tunately, from the relation between (A t ,t £ Z) and (At,t £ Z), it is easy to see, taking a multiplicative 
norm, that 7 (A) < j s (A), so that a sufficient condition for model (2.4) to possess a nonanticipative s.p.s. 
solution is that 7 s (A) < 0, which is the result of Theorem 2 of Bibi and Aknouche (2006). Nevertheless, 
using model (2.5) it seems to be difficult to show that the latter condition is also necessary. Using directly 
model (2.4), the following theorem shows that this condition is also necessary. 
Theorem 2.1 Model (2.4) admits a nonanticipative s.p.s. solution given by 



if and only if the top Lyapunov exponent 7 (A) given by (2.6) is strictly negative, where the above series 
converges a.s. for all t G Z. Moreover, the solution is unique and periodically ergodic. 

We now give a necessary condition for the strict stationarity, which coincide for the case S = 1 
(non periodic) to corollary 2.3 of Bougerol and Picard (19926) and which we will be needed in studying 
asymptotics for the QMLE. 

Corollary 2.2 // the P-GARCH model (2.1) possesses an s.p.s. solution then 



We now turn to an important consequence of strict periodic stationarity and the corresponding log- 
moment condition of Theorem 2.1, that is the existence of a moment of some positive order. A similar 
result for the non periodic GARCH case was proved by Nelson (1990) for the GARCH(1, 1) and Berkes 
et al. (2003) (Lemma 2.3) for the general GARCH case. 
Theorem 2.3 // 7 s (A) < then there is S > such that 




(2.8) 



k=l j=0 




where (3 t is the submatrix of A t defined by 





and 




(2.9) 
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3 Asymptotic properties of the quasi-maximum likelihood estimator 



Let 6 = (0 1 , 9 2 , . ., 9 S )' , where 9i = (o>j, a^i, a^q, f3 ii , /3 ip )', be the parameter vector which is supposed 
to belong to a compact parameter space C (]0,+oo[ x [0, +oo[( p+q ^ S . The true parameter value is 



unknown and is denoted by 9" = {9"{ ,9% , ...,9^)' with 9^ = (w?, a^, a£ g , $ tl , .... Consider a 

series {2/1,2/2; ■■■iVt-.T = NS} generated from model (2.1) with true parameter 9° 6 0, i.e. from the model 



2/f+Sn = \/K+S^Vv+Sn 

<? P 

VhSn = ^2 + E «° i^+Sn-i + E fijhv+Sn-j 
i=l j=l 



t G Z, 



(3.1) 



where the sample size T is supposed without loss of generality multiple of the period 5. The Gaussian 
logdikelihood function of 9 G G conditional on initial values yo, 2/i-<j, ho, h±- p is given (ignoring a 
constant) by 

1 N S 

L N s(e) = - w ^J2J2 l ^s(e), 



k=l s=l 



with 



Z t W = £L + log/ it , 
/it 



(3.2) 
(3.3) 



where ht is solution to the conditional model 



Vv+Sn 



hy+SnVv+Sn 
q 2 



h v+Sn =U>v + Yl OV2A,+S„-i + E (3 v ,jh v+ Sn-j 
i=l 3=1 



t > 1, 



(3.4) 



with given initial values yo, ...,yi- g ,ho, ...,hi- p . These values may be chosen, taking into account the 5- 
periodicity of the distribution of (yt), as an approximation of the non conditional variance. For instance 
for the first order P-GARCH model the non conditional variance is given by 



E(yt) =E(ht) 



S-lj-l 

ut + E EI («t-i + Pt-i) vt-j 
3=1 i=0 

5-1 

1 - 11 (<**-< + ^t-i) 

i=0 



which may be negative for some values of the strict periodic stationarity domain. Thus, taking at = (3 t = 0, 
we have the initial values 



Vq = h = uj , 
Vl-d=\-d = ^ 2 _ d , 



(3.5a) 
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which can be considered the same for the general P-GARCH model. Another choice of yo, y±- q , ho, 
hi-p is 

Vo = = Us, 

^=*-=«-ip (35i>) 

yi-d = hi-d = yfs-^ 

where d = max(p, q) and [k] = k if k > 1 and [&] = IS + k otherwise, I being the lower integer such that 
IS+k > 1. Obviously, these choices have only a practical value and do not affect the asymptotic properties. 
The QMLE of 9°, denoted by 9 ns, is the maximizer of Lns (9) on O and then the minimizer of 

1 N s „ 

fc=l s=l 

Let 7 5 (A ) denote the top Lyapunov exponent associated with (A^t G Z) where yl^ is just A t defined in 
(2.4) with 9° in place of 9. To study the strong consistency of 9ns consider the following assumptions. 

Al: 7 5 (.4°) < and V# G 9, p ( f[ /O < 1. 

i p 

A2: The polynomials a® (z) = Yl ct^jZ^ and j3 v (z) = 1 — Yl P^jZ^ have no common roots, a° (1) / 

j=i ' i=i 
and a° q + /?° p / for all 1 < v < S. 

A3: (r] t ) is non degenerate. 

As seen in Theorem 2.3, assumption Al ensures the existence of a finite moment, for the solution of 
(3.1), which is the key for proving consistency and asymptotic normality irrespective of any moment condi- 
tion. The condition p ^ Y\ Pv*j < 1 is imposed for any 9 G O in order to obtain h t (9) as a causal solution 
of {yt,yt-i, •••} while A2 is made to guarantee the identifiability of the parameter. Similar conditions for 
the GARCH case have been considered by Francq and Zakoi'an (2004). 

Given a realization {yi, ...,yNs}, h (#) can be approximated for 1 < t < NS by 

k(9) = j^ f) +logh t (9), (3.6) 

where ht (9) is the s.p.s. solution of 

q p 

K+Sn = V v + X Otv,iyl+Sn-i + X Pv,jh v +Sn~j, t G Z, (3.7) 
i=l j=l 

for 9 G with h t = h t {9°). Let 

N S 

L N s{e) = -—Y.Y, l ^{0). 

k=l s=l 
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The following theorem establishes the strong consistency of 9 ns- 
Theorem 3.1 

Under A1-A3 9ns is strongly consistent in the sense that 9ns — * 8° a- s - as N — > oo. 
Remarks 

1) From assumptions A1-A3, it is clear that the above result remains true for the particular periodic 
ARCH (P-ARCH) case, i.e. when (3 tJ = 0. 

2) In order that Al holds for, for example, the particular P-GARCH(1, 1), we may chose the parameter 
space as a compact of the form = ([e, 1/e] x [0, 1/e] x [0, 1 — e]) S where e > is so small that the 

parameter 9° = (oj®, a\ 1 , /3§ 1; u®, 1 , &\ ]_..., Wg, ajg v /3g x )' belongs to 6. The fact that /?° 1 must be lower 

S 

than 1 for all v is so restrictive and can be weakened so that f\ P%,i < 1 which is just assumption Al for 
the particular P-GARCH(1, 1). 

Now we turn to the asymptotic normality of 9 ns- Consider the following additional assumptions. 
A4: 9° is in the interior of 0. 
A5: E {rft) < oo. 

Assumption A4 is standard and allows to validate the first order condition on the maximizer of the 
log-likelihood while A5 is necessary for the existence of the limiting covariance matrix of the QMLE. 

The following result establishes the asymptotic normality of 9 ns- 
Theorem 3.2 
Under A1-A5 we have 

VNS (9 NS -6°)~*N (0, (E {vt) - 1) J~ l ) , as TV - oo, (3.9) 
where the matrix J given by 



J :-- 



v=l 



3989' 



v , 1 8h v (9°) d_K (0°) i , . 

^ Eec \W¥)~de (3 - 10) 



v=l 



is block diagonal. 

Let us now apply the forgoing results to the first order P-ARCH process given by 



Vv+Sn = Vv+Sn\l^v + a vVv+Sn-V 

where u;° > and a° > 0. It is easily seen that the strict periodic stationarity condition for the P-ARCH{\) 
reduces to 



< J] a° v < exp (-E (log ( V 2 V ))) := a, 

v=l 

under which, supposing that 9° = (m®, a\ 1 ,oj^ a 2 u ■■■■> ^s^ a s l)' belongs to a compact of the form 

S 

= ([6,1/e] x [0y/ 5 -l]) , (3.11) 
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the QMLE is then by Theorem 3.1 strongly consistent. 

Moreover, if we assume that 9° is in the interior of given by (3.11), then from Theorem 3.2 the 
QMLE is also asymptotically Gaussian and its limiting distribution is given by (3.9) and (3.10) which in 
the P-GARCH(1, 1) case reduces to 

/ Jl 02x2 • • • 02x2 \ 

02x2 J2 ■ ■ ■ '■ 

■ '■ '• 02x2 

\ 02x2 • • • 02x2 Js ) 



J 



where 



J v — E 



, 1 < v < S. 



4 Appendix. Proofs 

4.1 Appendix A. Proofs of Section 2's results 
4.1.1 Proof of Theorem 2.1 

We only give a proof of the necessity part of the theorem as the sufficiency one was established earlier by 
Bibi and Aknouche (2006). Suppose that model (2.4) admits a nonanticipative s.p.s. solution (Y t ,t G Z). 
Then iterating (2.4) m -times for some m > 1 and some v G {1, S} gives 

m j — 1 m+1 

Y v = ^ ^ J J 4,- i^'r j ~t~ • b- j ^ V I • 

j=0 i=0 j=0 

and exploiting the non-negativity of the coefficients of At and those of Y v it follows that for all m > 1 

m j-1 

Y v > ^^ || A v —jB v —ji a.s. 

J- 1 



implying that the series \ , Y\ A v -iB v -j converge a.s. and therefore 



j=0 



3-1 



|~J A v _iB v _j — > 0, a.s. as j — > 00, 



from which we want to show that 



j'-i 



I J Aj_j — > 0, a.s. as j — > 00. 



(Al) 



(A2) 



i=0 
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This holds whenever 



lim TT A v -iek = 0, a.s. for all k G {1, r} , (^4-3) 

J i=0 

where (e^) is the canonical basis of W. Observing that B v -j = io v -jf^ v _^e\ +oj v -je q+ i, (A.l) implies that 

j'-i j'-i 
^v-jvt-j Y\_ A v -i&i — ► 0, a.s. as j — > oo, and u;„_j A v _ie q+ i — ► 0, a.s. as j — > oo, (^4-4) 
i=0 i=0 

and since o;„_j > 0, this shows (A. 3) for k = q + 1. Moreover, since for all A; = 1, ...,p 

A v ~j+ie q+ k = P v _j +lk rj v _j +1 ei + v _j +lk e q+ i + e q+ k + i, 
then for k = 1 we have almost surely 

= lim TT^_je ?+ i > lim TT A v -ie q+2 > 0, 

?'— >oo -*- j— >oo -*■ 

i=0 i=0 

showing (A. 3) for k = q + 2 and then recursively for fc = q + j, j = l...p. On the other hand, since 

A v - j+ ie q = a v -j+i jq ril_j +1 ei + a v - j+1:q e q+1 , 
then from (A4), follows (A3) for k = q. Finally, from the identity 

A v -j + \ek = tt v ~j+i,kvl-j+i e i + a v-j+i,k e q+i + e fc+i) fc = 1, ...,q — 1, 
follows (A3) for the rest values of A; using a backward recursion. This show (A2), which implies that any 
subsequence of the sequence Uj := Y[ A-« converge a.s. to zero as j — > oo. In particular, 



«=o 



Sj-i j-i /S-i \ 

C/ 5i = JJ = JJ I Y[ Aj-iS-k I -> 0, a.s. as j -> oo. (A5) 

i=0 i=0 \k=0 J 

fS-l \ 

Now, since the sequence I Q Aj-iS-fc>2 G Z J is i.i.d. for all i>, then from (A.5) and Lemma 2.1 of 

Vfc=o / 

Bougerol and Picard (19926) it follows that the top Lyapunov exponent associated with the i.i.d. sequence 

5-1 \ 

f7 A v -i,s-k,i G Z I, which is exactly 7 5 (A), is strictly negative. This completes the proof. 
k=o J 

4.1.2 Proof of Corollary 2.2 

As the (At) is nonnegative, the corresponding top Lyapunov exponent 7 5 (A) is greater that the top 
Lyapunov exponent corresponding to nonrandom sequence of matrices, say (3 = (f3 t ,t G Z) obtained when 
one set the elements of the first q lines and those of the first q columns to be equal to zero. In other words, 



7 5 04) >7 5 (/3) :=logp^n^ 
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Thus, if the model has an s.p.s. solution which is equivalent to the condition 7 s (A) < 0, the conclusion of 
the Corollary is true. 

4.1.3 Proof of Theorem 2.3 

The proof is similar to that of Lemma 2.3 of Berkes et al. (2003). First we have to show that if 7 s (A) < 
then there is 5 > and no such that 



E (\\A nQS A no s-i...A 1 \\ s ^ < 1. 



(A6) 



Since j s (A) = inf ne pj. {^E (log ||j4 T jS^4„s_i...Al||)} is strictly negative, there is a positive integer no such 
that 

E (log\\A noS A noS -i...A 1 \\) < 0. 

On the other hand, working with a multiplicative norm and by the property i.p.d. of the sequence (A t ) we 
have 



E (\\A noS A noS _ 1 ...A 1 \\) = \\E(A noS A noS ^ 1 ...A 1 )\\ 

= \\E(A s A S - 1 ...A 1 ) no \\ 
< WEiAsAs-L-.A^r <oo. 

Let f(t) = E(\\A noS A noS -i...A 1 \\ t ). Since /' (0) = E (log || A noS A noS -i...A 1 1|) < 0, / (*) decrease in a 
neighborhood of and since / (0) = 1, it follows that there exists < 5 < 1 such that (A.6) holds. 
Now from Theorem 2.1 we have for some v £ {1, S} 



\Y V \\<T, 



k=i 



and since < 5 < 1 then 



k=i 



k-l 
j=0 



k-i 

3=0 



I E v — ^|| -|- Iji^^ll , 



|-E>D-fc|| 5 + \\B V \\ S , 



which by the independence of A v -j and B v _k for j < k implies 



00 


( 


k~i 


5 \ 


E\\Y V \\ S < J2 E 








k=i 


V 


j=0 


) 



E (\\B v ^ k \\ S ^ + E (\\B v \\ 5 ^j 



< B(5)J2E 



k=i 



k-l 
j=0 



+ E 



(\\B V \\ 5 ), 



12 



where B (6) = max E (\\B v ^k\\ 5 )- Using (A.6) there exist a v > and < b v < 1 such that 

o<v<s-i V / 





( 


k-l 




E 






<a v b k < ab k 




{ 


3=0 





where ab k = max {a v b k \. This prove that E \\Y V \\ 5 < oo and then (2.9). 

0<V<S~1 

4.2 Appendix B. Proofs of Theorems 3.1-3.2 

The proofs of Theorems 3.1 and 3.2 are by now standard and follow from similar arguments used in showing 
the strong consistency and asymptotic normality of the QMLE for standard GARCH models (Francq and 
Zakoi'an, 2004). The main aim here is to reveal the basic assumptions and to quantify the asymptotic 
distribution of the QMLE for the P-GARCH . Since there are several similarities between the standard 
GARCH and the P-GARCH , certain steps of the proof for the P-GARCH case are similar in spirit to 
that of the standard GARCH one. Thus, we give details of proof only when it seems pertinent and refer 
to Francq and Zakoi'an (2004) for the details. 



4.2.1 Proof of Theorem 3.1 

Theorem 3.1 will be proved by showing several lemmas. Lemma B.l establishes the uniform asymptotic 
forgetting of initial values, Lemma B.2 ensures identifiability of the parameter, Lemma B.3 shows the finite- 
ness of the limiting criterion J2s=i Ee° (#)) an d that this one is uniquely minimized at 9°, while Lemma 
B.4 uses the compactness of 9 and a periodically ergodic argument to conclude the strong convergence. 
Lemma B.l Under Al we have 



sup 

eee 



L NS (6) - L NS (9) 
Proof Rewriting (3.7) in a vectorial form as follows 



0, a.s. as N — > oo, 



(A.7) 



where h t = (h t , h t -i, h t - v +i)' and a t = (oj v + J2 a l) , i ^ +5n _ i , 0, j . From Corollary 2.2, the 

V i=l / lxp 

assumption p (^\\ (3^j < 1 of Al and the compactness of we have 



supp ( TT/3« 1 < i, 



(A.8) 
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which implies, by iterating (A 7), that 

t-i fe-i 



k=0 i=0 



i=0 



If we denote by /i t and « t the vectors obtained from h t and a t , respectively, while replacing h t -j by ht-j 
with initial values given by (3.5), then we have 



t-i fc-i 



t-g-1 fc-l 

k=0 i=0 k=t-q i=0 i=0 



and from (A8) it follows that 



sup 

eee 



ht-ht 



= sup 

eee 

< Mp\ 



t-l k-l 



k=t—q i=0 



i=0 



from which, we get 

L NS (0) - L NS (9) 



NS 



sup 



1 



t=i 



eee 



h t - h t 



h t h t 



y 2 t + 



log 



< 



max sup to 

l<v<S eee 



where the inequality llog a | < t for positive x and u has been used. The existence of a moment for 

yt (cf, Corollary 2.2) implies, by the Borel-Cantelli lemma, that p l yl — > a.s. so that the conclusion of the 
lemma follows from the Toeplitz lemma. 

Lemma B.2 Under A1-A3 there is t G Z sitc/i tftai ?/ fo t (0) = h t (9°) a.s. then 9 = 9°. 

Proof From the assumption p ^ f\ Pvj < 1 of Al the polynomials (/3° (L)) 1<v<s (L being the backward 

shift operator) are invertible. Suppose that h t (9) = h t (9°) a.s. for some t, then using (3.7) we have 



oty (L) _ ag (Ly \ 2 
Pv(L) (3° V (L) 



Vv+St 



for all 1< v < S. 



some 1 < v < S then there exists a deterministic periodic time-varying combination 

Pv ( L ) Pv(L) ~ ~ 

of Vv+st~j> 3 ^ L This contradict the fact that y 2 v+st = E {vl+st/vl+st-n •••) + VhSt (^+St ~ *)> since 
by A3 (r) t ,t € Z) is non degenerate. Therefore, 

a v (z) a° v (z) . . w° w„ 

- " V \z \ < 1 and — ^ . - , for all 1 < v < 6, 
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which by the assumption A2 of absence of common roots implies that a v (z) = a° (z), (3 V (z) = /?° (z) and 
lo v = uj® for all 1 < v < S, proving the lemma. 
Lemma B.3 Under Al 

s 

J2 E <>° ( l * f))<^ 

v=l 

and Ylv=i E s° (lv {&)) is minimum at = 9°. 
Proof By Theorem 2.3 we have 



s s 
Y,E e o(\ogh v (9 )) = ^TEe* 1 - (log h v (9°) S ) 



v=l 



v=l 



1 s 

< ^logEeo (h v (9°) S ) <oo, 
ii=i 



from which it follows that 



v=l 



v=l 



hv (0°) V V 



i^ + iogh v f) 



K f) 

S + J2 E 0° ( l °Sh v (0°)) < oo. 



v=l 



Finally, 



v=l 



v=l 



Y^Eeo{i v (e))-Y,E e o{i v {e°)) = ^TE 6 o 

v=l 

s 



v=l 



, K{9) Kf) 
g K(e°) + K{6) 

h v (g) , K (g°) 

K \ f) K (9) 



0, 



showing that the limit criterion is minimized at 9°. 

Lemma B.4 Under Al, for all 9^9° there is a neighborhood V (9) such that 



(A9) 



lim inf _ inf ( ——Lns 

n^oo eev(e) \ N 



> 



v=l 



Proof For all 9 £ and all integer k, let Vk (9) an open sphere of center 9 and radius 1/k. Using Lemma 
B.l we have 



lim inf_ inf 



limnj^no inf inf;; 



{-JjLns (*)) 



1~ /~\ \ lim A r ^oo lilil±li e e Vfe(e)ne 

limjv-oo sup jj sup eee L NS {9) - L NS (9) 

1 JV_1 5 „ 
> lim inf— VV, inf l v+nS ((?) . 
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Applying the ergodic theorem for the i.i.d. sequence ^Ec=i^+nS^jj with E ^X)f=i l v +nS (j^j ) £ 
MU {00} (cf, Billingsley 1995, p. 284, 495) it follows that 

liminf— y^y] inf l v+nS ( 9 ) = V] inf Z^+nS ( 

and by the Beppo-Levi theorem (e.g. Billingsley, 1995, p. 219), we have 

S ( \ 5 

y2 E 6° ~ illf l v+nS ( O) I Yl E 0° W) as /C ^ CX), 

^ Veev fe (e)ne wy ^ 
which by (A 9) proves the lemma. 

Proof of Theorem 3.1 In view of Lemmas B.1-B.4, the proof of the theorem is completed by using an 
argument of compactness of 0. First, for all neighborhood V (0°) of 0° we have 



lim sup inf ( -1~L NS fa)) < lim ( -^-L NS (9°)) = lim (-^- T L NS (9°) 

s 

v=l 

The compact is recovered by a union of a neighborhood V (0°) of 9° and the set of neighborhoods V (9), 
9 £ 0\V (9°), where V (0) fulfills Lemma B.4. Therefore, there exists a finite sub-covering of by V (0o) , 
V(0i),...,V(0 fc ) such that 

inf ( ——Lns (9)) = min inf ( — — Ljvs M • 

From (A. 10) and Lemma B.4, the latter relation shows that 9ns G V (0°) for A sufficiently large, which 
complete the proof of the theorem. 

4.2.2 Proof of Theorem 3.2 

The proof of Theorem 3.2 is based on a Taylor expansion of ~ at 9°, along the lines of Francq and 



Zakoian (2004), which is given by 

NS rCT NS / NS r, 2 l 

° = E § M = W" 1/2 E § «n + W E 5^ («) ("I 1 ' 2 (fc™ - «°) . 



where the coordinates of are between the corresponding entrees of ns and those of 0° . 

The proof is also based on several lemmas, Lemma B.5-B.9, which are aimed at establishing: the 
integrability of the first derivatives of the limiting criterion at 0°, the invertibility of J and its relation 
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with de first derivatives of the limiting criterion, the uniform integrability of the third-order derivatives of 
the limiting criterion, the asymptotic forgetting of starting values for the derivatives, and a central limit 
theorem for martingales differences together with a periodically ergodic theorem for the second derivatives 
of the criterion. 
Lemma B.5 We have 

S r\i //in\ Ai / rtr>\ $ 



s=l 



dl s (0°) dl s (6°) 



89 89' 



< 



oo, and Ego 



s=l 



8% (9°) 



8989' 



< oo. 



Lemma B.6 Under A1-A5 J is invertible and 

S „, ,„„x 2 



s=l 



dl, 



89 



(E (r,t) ~ 1) J- 



Lemma B.7 The following limit relations 



N S 



k=l s=l 



dh+ks (9°) 8l s +ks (9°) 



89 



and 



N S 

sup N-^s^y] 

Oev{0°) ' * ' * 



k=l s=l 



81 



s+kS 1 



89 



dls+ks (9°) 



o P (1) , 



89 



89 



o P (1) , 



are true. 

Lemma B.8 There is a neighborhood V (0°) of 9° such that for all i,j, k E {1, S(p + q + 1)} 



E^o sup 
,=i eev(e") 

Lemma B.9 The following limit results are true 



8\ (9) 



89 i d9 j d9 k 



< oo. 



N S 



^EE^ilr * N (°- ^ W) " x ) J ) and £ £ w ( ? ) ^ J < 



k=l s=l 



fe=l s=l 



The proofs are very similar to those of Francq and Zakoi'an (2004). It suffices to replace the stationarity 
and ergodicity arguments by the periodic stationarity and periodic ergodicity ones, respectively. For this 
we omit the proofs which are quite lengthy. 
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